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ABSTRACT 

We study thermal structure and evolution of magnetars as cooling neutron stars with 
a phenomenological heat source in a spherical internal layer. We explore the location 
of this layer as well as the heating rate that could explain high observable thermal 
luminosities of magnetars and would be consistent with the energy budget of neutron 
stars. We conclude that the heat source should be located in an outer magnetar's 
crust, at densities p<5x 10 11 g cm -3 , and should have the heat intensity of ~ 10 20 
erg cm~ 3 s . Otherwise the heat energy is mainly emitted by neutrinos and cannot 
warm up the surface. 

Key words: dense matter — stars: magnetic fields — stars: neutron - neutrinos. 



1 INTRODUCTION 

There is growing evidence that soft gamma repeaters 
(SGRs) and anomalous X-ray pulsars (AXPs) belong to 
the same class of objects, magnetars, which are warm, iso- 
lated slowly rotating neutron stars of age t < 10 5 yr with 
unusually strong magneti c fields, B > 10 14 G (see, e.g., 
IWoods fc ThompsonlEool. for a recent review) . There have 
been attempts to explain the activity of these sources and 
the high level of their X-ray emission by the release of the 
magnetic energy in their interiors but a reliable theory is 
still absent. 

In this paper we analyze the thermal evolution of mag- 
netars as cooling isolated neutron stars. We do not attempt 
to construct a self-consistent theory of the magnetars but 
instead address the problem phenomenologically. We show 
that magnetars are too hot to be treated as purely cooling 
neutron stars; they require some heating source, which we 
assume operates in their interiors. Our aim is to analyze the 
location and power of the heating source that are consistent 
with observed thermal luminosities of SGRs and AXPs and 
with the energy budget of an isolated neutron star. 



2 PHYSICS INPUT 

We use our general relativistic cooling code (Gnedin, 
Yakovlev & Potekhin 2001). It simulates the thermal evo- 
lution of an initially hot isolated neutron star taking into 
account heat outflow via neutrino emission from the star 
and via thermal conduction within the star, with subse- 



quent thermal photon emission from the surface. To fa- 
cilitate calculations, the star is artific ially divided (e.g., 
iGudmundsson. Pethick fc Epstein! Il983l) into a thin outer 
heat blanketing envelope (extending from the surface to the 
layer of density p = pb ~ 10 — 10 11 g cm -3 , with the thick- 
ness of a few hundred meters) , and the bulk interior (from pb 
to the center). In the bulk interior the code solves the full set 
of equations of thermal evolution in the spherically symmet- 
ric approximation, neglecting the effects of magnetic fields 
on thermal conduction and neutrino emission. In the blan- 
keting envelope the code uses the solution of the stationary 
thermal conduction problem obtained in the approximation 
of a thin plane-parallel layer for a dipole magnetic field con- 
figuration. This solution relates temperature Tb at the base 
of the blanketing envelope (p = pb) to the effective surface 
temperat ure T s properly averaged over the neutron star sur - 
face (e.g.. |Potekhin fc Yakovlevl200ll.|potekhin et ai]|2003l) . 



In the code we have mainly used the Tb — T s rela- 
tion obtained specifically for the present study, assuming 
pb = 10 10 g cm -3 and the magnetized blanketing envelope 
made of iron. Magnetars are hot and have large tempera- 
ture gradients extending deeply into the heat blanketing en- 
velope. Thus even high magnetic fields do not dramatically 
affect the overall thermal conduction in th e envelope and 
the T b — T B relation (as discussed, e.g., bv lPotekhin et all 
2003). Since magnetars are hot inside, these fields should 
not greatly affect heat transport in the bulk of the star (at 
p > pb). This justifies the approximation of spherically sym- 
metric temperature distribution at p > pb- At lower temper- 
atures the electron thermal conductivity in the interior be- 
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comes strongly anisotropic and such an approximati on can 
be questionable llGeppert. Kiiker fc Pagdl2004 l200dL 

Our standard cooling code includes the effects of mag- 
netic fields only in the blanketing envelope. In this re- 
spect our present model of the blanketing envelope with 
p b = 10 10 g cm -3 seems less adeq uate than the previou s 
model with p b = 4 x I0 11 g cm" 3 JPotekhin et al]l2003lh 
However, the latter model, by construction, neglects neu- 
trino emission from the blanketing envelope. Because the 
neutrino emission from the layers with p > 3 x 10 10 g cm" 3 
in hot magnetars is important (Section |3J, the blanketing 
envelope with pb = 10 10 g cm -3 is expected to be more ap- 
propriate for the present investigation. That is why we have 
performed calculations with this latter model but have done 
additional tests using the model with pb = 4 x 10 11 g cm -3 
(Section |OJ. 

We have also included, in a phenomenological manner, 
the effects of magnetic fields on the thermal evolution of the 
bulk stellar interior. Most importantly, we have introduced 
a heat source located within a spherical layer in the interior 
(it can be associated with the magnetic field, Section [4.211 . 
The heating rate H [erg cm -3 s _1 ] has been taken in the 
form 



H = H Q{p u p 2 ) exp(-t/r) 



(1) 



where Hq is the maximum heat intensity, 0(pi, p2) is a step- 
like function (O ~ 1 within the density interval pi < p < p 2 \ 
B ~ outside this interval, with a sharp but continuous 
transitions at the interval boundaries), t is the star's age, 
and t is the life time of the heating source. A specific form 
of H is not important for our main conclusions. We do not 
specify the nature of this source. We set t = 5x 10 4 years to 
explain high thermal states of all observed SGRs and AXPs 
(Section|3Jl. We do not consider longer r which would require 
higher energy budget (while the budget is already severely 
restricted even for r = 5 x 10 4 years). We treat Ho, pi and 
p2 as free parameters with the aim to understand what the 
intensity and the location of the heat source should be in 
order to be consistent with observations and with the energy 
budget of an isolated neutron star. 

It is instructive to introduce the total heat power W 00 
[erg s _1 ], redshifted for a distant observer, 



W°°(t) = J dVt 



(2) 



where dV is the proper volume element, and <& is the redshift 
metric function. 

In the neutron star core we use the equa- 

tion of state of dense mat ter constructed by 

lAkmal. Pandharipande fc Ravenhalll (|T998) (model Argon 
V18 + 8v + UIX*); it is currently considered the most elab- 
orated equation of state of neutron star matter. Specifically, 
we employ a conve nient parameterization of this equat ion of 
state proposed bv iHeiselbe^g ^j^Frjo^th- Jensenl (1999) and 
described as APR III bv lG^sako^e^LT(l200^) . According 
to this equation of state, neutron star cores consist of 
neutrons, protons, electrons, and muons. The maximum 
(gravitational) neutron star mass is M = 1.929 Mq. 
The powerful direct Urea process of neutrino emission 
llLattimer et al]ll99ll) is allowed only in the central kernels 
of massive neutron stars with M > 1.685 Mq (at densities 
p > 1.275 x 10 15 g cm" 3 ). 
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Figure 1. Observational data on the blackbody surface tem- 
peratures T s °° of seven magnetars. The shaded rectangle is the 
"magnetar box". The data are compared to the theoretical cool- 
ing curves of the 1.4 Mq neutron star with B = 5 x 10 14 G and 
no internal heating, either without superfluidity (the solid line) 
or with strong proton superfluidity in the core (the dashed line 
SF). 



We use neutron star models with two masses, M = 
1.4 M e and M = 1.9 Mq. The 1.4 Mq model is an exam- 
ple of a star with standard (not too strong) neutrino emis- 
sion in the core (provided by the modified Urea process 
in a non-superfluid star). In this case the (circumferential) 
stellar radius is R — 12.27 km, and the central density is 
p c = 9.280 xlO 14 gem" 3 . The 1.9 Mq model (R = 10.95 km, 
p c = 2.050 x 10 15 g cm -3 ) is an example of a star whose neu- 
trino emission is enhanced by the direct Urea process in the 
inner core. 

We have updated the thermal conductivity of electrons 
and muons in the stellar core by new results (P. S. Shternin 
& D. G. Yakovlev 2006, in preparation) which take into ac- 
count the Landau damping of transverse plasmons in the in- 
teraction of electrons and muons with surrounding charged 
particles (following the resu lts for quark plasma obtained by 
Hcisclbcrg fc Pethickll993l) . This update has not noticeably 
affected our cooling scenarios. 



3 RESULTS 

For the observational basis we take seven sources (two 
SGRs and five AXPs indicated in Figure 0. The esti- 
mates of their spindown ages t, surface magnetic fields 
B and the blackbody surface temperatures (red- 
shifted for a distant observer) are taken from T ables 
14.1 a nd 14.2 of the review pa per bv IWoods fc Tho mpson! 
12006ft and from the paper bv lMcGarrv et al] J2005f). The 



data are from the origi nal publicatio n s of iKulkar ni et alJ 
feOOSt CS GR 0526-66): IWoods et alJ j200ll. 120021) fSGR 



1900+14); iGotthelf et alJ feOOSft and 



Morii et al 
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Table 1. Four positions of the heating layer, and the heating power W°° for the 1.4 Mq star with Ho = 3 X 10 20 erg cm 3 s 1 and 
t = 1 kyr 



No. pi (g cm 3 ) p2 (g cm 3 ) W°° (erg s 1 ) 

I 3 X 10 10 10 11 4.0 X 10 37 

II 10 12 3 x 10 12 1.9 x 10 37 

III 3 x 10 13 10 14 1.1 x 10 38 

IV 3 x 10 13 9 x 10 14 1.1 x 10 39 



(IE 1841 -0451: iMcGarrv et all feOOff) (CXOU J010043.1- 
721134); IGavriil fc KaspiT 2002ll and iRea. et al.l EoB 
(1RXS _J1 70849^-400910); IGavriil fc Kaspil j2002ll and 



IPatel et, alJ (120031) (4U 0142+611: IGavriil fc KaspTptua ) 
and IWoods et alJ l)2004h (IE 2259+586). In the absence of 
T a °° estimates for SGR 1627-41 we do not include this SGR 
in our analysis. We have also excluded SGR 1806-20 and 
several AXPs whose thermal emission component and char- 
acteristic age seem less certain. The radiation from the se- 
lected sources has the pulsed fraction < 20%. This indicates 
that the thermal radiation can be emitted from a substantial 
part of the surface. 

The blackbody surface temperatures of the selected 
sources are plotted i n Fig ure versus spindown ages t. 
IWoods fc Thompson! l)2006h present these data without for- 
mal errors, which are actually large. We introduce, some- 
what arbitrarily, the 30% uncertainties into the values of 
r a °° and the uncertainties by a factor of 2 into the values 
of t. The data are too uncertain and our cooling models 
are too simplified to explain every source by its own cooling 
model. Instead, we try to explain the existence of magne- 
tars as cooling neutron stars that belong to the "magnetar 
box", the shaded rectangle in Figure (attributed to an av- 
erage persistent thermal emission from magnetars, excluding 
bursting states). Our results will be sufficiently insensitive 
to the neutron star mass, and we will mainly use the neu- 
tron star model with M = 1.4 Mq (unless the contrary is 
indicated). Similarly, the results will be not too sensitive to 
superfiuid state of stellar interiors and we will mostly neglect 
the effects of superfluidity of nucleons in the stellar crust and 
core. To be specific, we mainly assume the dipole magnetic 
field in the blanketing envelope with B = 5 x 10 14 G at 
the magnetic poles. Some variations of B will not affect our 
principal conclusions (Section 14.11 . 

In Figure Q we show the theoretical cooling curves 
T a °°(i) for the 1.4 Mq isolated magnetized neutron star 
without any internal heating. The solid line refers to a 
nonsuperfiuid neutron star while the dashed line SF is 
for strong proton superfluidity in the stellar core. This 
superfluidity strongly suppresses neutrino emission in the 
core and n oticeably increases r „ °° at the neutrino cooling 
stage (e.g.. lYakovlev fc Pethickl 120041) . Let us stress that 
the surface temperature of these stars is highly nonuni- 
form, with the magnetic poles being much hotter than the 
equator. In th e figures we plot the average surface tem- 
perature (e.g., IPotekhin et al-llioOSt) . Clearly, the magne- 
tars are much hotter than ordinary cooling neutron stars. 
The observations of ordinary neutron stars can be ex- 
plained by the cooling theory without any reheating (e.g., 



lYakovlev fc Pethickl 120041: iPaee. Geppert fc Weberl [2006), 
while the observations of magnetars suggest that the mag- 
netars have additional heat sources. We assume further that 
these sources are located inside magnetars. According to al- 
ternative models, powerful energy sources can be available in 
magnetar magnetospheres l Tjiompson fc Beloborodovl2005l : 



Bclo borodov fc Thompson! 20061. 



We introduce the internal sources in a phenomenologi- 
cal way described in Section All results presented below 
(Figures are obtained including the internal heating in 
accordance with Eq. Q. They indicate that the magnetars 
are hot inside, with the temperature T ~ 10 9 K (or even 
higher) in the crust, at p > 3 x 10 10 g cm" 3 . Such stars are 
very strong sources of neutrino emission, which is vitally im- 
portant for the magnetar thermal structure and evolution. 

Our simulations of cooling neutron stars with a pow- 
erful internal heating show that after a short initial relax- 
ation (t < 10 years) the star reaches a quasi-stationary state, 
which is fully determined by the heating source. The energy 
is mainly carried away by neutrinos, but some fraction is 
transported by thermal conduction to the surface and radi- 
ated away by photons. The interior of these cooling neutron 
stars is highly non-isothermal. The hotter layers are those 
where the heat is released and where the neutrino emission 
is not too strong; these layers are located in the crust. The 
approximation of isothermal interior, which is excellent for 
ordinary middle-aged neutron stars, cannot be used while 
studying the cooling of magnetars. 

The left panel of Figure |5] shows the temperature dis- 
tribution inside the 1.4 Mq star of age t — 1000 years. This 
distribution remains the same during the entire magnetar 
stage (t < r, see Eq. Q) and we have chosen t — 1000 years 
just as an example. We have considered four locations of the 
heat layer, p\ — p2, summarized in Tabled (I) 3 x 10 10 — 10 11 
g cm -3 (in the outer crust, just below the heat blanketing 
envelope), (II) 10 12 - 3 x 10 12 g cm" 3 (at the top of the 
inner crust), (III) 3 x 10 13 - 10 14 g cm" 3 (at the bottom of 
the inner crust), and (IV) 3 x 10 13 - 9 x 10 14 g cm" 3 (at the 
bottom of the inner crust and in the entire core). These loca- 
tions are marked by hatched rectangles. The density ranges, 
which are appropriate to the outer crust, the inner crust and 
the core, are indicated in the upper part of the figure. Let 
us remind the reader that the outer crust has a thickness of 
a few hundred meters and a mass of ~ 1O" 5 M , the inner 
crust can be as thick as 1 km and its mass is ~ 10" 2 Mq, 
while the core is large (radius ~ 10 km) and contains ~ 99% 
of the stellar mass. Therefore, the heating layers I, II, and III 
are relatively thin, while the heating layer IV is wide and 
includes most of the stellar volume. We have allowed our 
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Figure 2. Left: Temperature profiles within the 1.4 Mq neutron star of age t = 1000 years with four different positions I— IV of the 
heating layer (given in Table IT1 and indicated by hatched rectangles) and two levels of the heat intensity Ho = 3 X 10 19 and 3 X 10 20 
erg cm~ 3 s _1 . The magnetic field is B = 5 X 10 14 G. Right: Cooling curves for these models compared with the observations. 



heat sources to have two different intensities, Ho = 3 x 10 
and 3 x 10 20 erg cm -3 s _1 . For illustration, in Tabled we 
present also the heating power W°° calculated from Eq. 
for the four layers in the 1.4 Mq star of age t = 1000 years 
at H = 3 x 10 20 erg s _1 . 

In all the cases the neutron star core appears to be much 
colder than the crust because the core quickly cools down 
via strong neutrino emission (via the modified Urea process 
for the conditions in Figure |5J. Placing the heat sources far 
from the heat blanketing envelope is an inefficient way to 
maintain warm surface; the heating layer can be hot, but 
the energy is radiated away by neutrinos and does not flow 
to the surface. For the deep heating layers (cases II, III, or 
IV), the heat intensity Ho = 3 x 10 19 erg cm -3 s _1 is clearly 
insufficient to warm the surface to the magnetar level. The 
higher intensity 3 x 10 20 erg cm~ 3 s _1 helps but is still less 
efficient in these layers than in the layer I, which is close to 
the heat blanketing envelope. For the latter Ho the heating 
of the crust bottom (case III) and the heating of the entire 
core (case IV) lead to the same surface temperature of the 
star. Thus, the most efficient way to warm the surface is 
to place the heating layer in the outer crust, near the heat 
blanketing envelope. This conclusion is further supported by 
Figures 

The right panel of Figure^ shows cooling curves of the 
1.4 Mq stars for the same models of the heating layer as 
in the left panel (with one exception — we do not show 
the cooling curves for the case IV, for simplicity). The cool- 
ing curves of the star of age t < 5 x 10 4 years are almost 
horizontal, indicating that these stars maintain their high 
thermal state owing to the internal heating. Non-horizontal 
initial parts (t < 100 years) of two curves, which correspond 
to deep and non-intense heating, show the residual initial 
relaxation to quasi-stationary thermal states. One can see 
that the only way to explain the sources from the "magnetar 
box" is to place the source in the outer crust and assume 




1 2 3 4 5 6 
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Figure 3. Same as in the right panel of Figure El for the three 
positions of the heating layer, one value of the heat intensity and 
two neutron star masses, 1.4 and 1.9 Mq. 

-ffo ~ 10 20 erg cm~ 3 s" 1 . As the heating is exponentially 
switched off at t > 5 x 10 4 years, following Eq. l|T). the 
surface temperature drops down accordingly. The magnetar 
stage is over and the star transforms in to the ordinary cool- 
ing neutron star jPotekhin et al.ll2003T) which cools mainly 
via the surface photon emission. 

Figure compares cooling curves of neutron stars of 
two masses, 1.4 Mq and 1.9 Mq, for the same three loca- 
tions I— III of the heating layer and for one heat intensity 
H = 3 x 10 20 erg cm -3 s -1 . The thick lines are for the 
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1.4 Mq star; they are the same as in the right panel of Fig- 
ure|S] The thin lines are the respective curves for the 1.9 Mq 
star, which is a very efficient neutrino emitter because of 
the direct Urea process in its core (Section |5J. If the heat- 
ing layer is in the location II or III (in the inner crust), the 
enhanced neutrino emission from the core of the massive 
star carries away a substantial amount of heat and decreases 
the surface temperature of the star. If, however, the heat- 
ing layer is placed in the outer crust (case I), direct Urea 
in the core has almost no effect on the surface temperature. 
In this case the surface temperature is almost insensitive to 
the physical conditions in the stellar core and in the inner 
crust, in particular, to the neutrino emission mechanisms 
and superfluid state of matter. In other words, surface lay- 
ers are thermally decoupled from the deep interior. A similar 
situation occurs in ordinary young and hot cooling stars in 
the initial cooling stage, before internal thermal relaxatio n 
jLattimer et aUll994t iGnedin. Yakovlev fc PotekhinllioOll) . 
This justifies our neglect of the effects of superfluidity in 
the calculations (although these effects can be vitally im- 
portant for ordinary middle-aged cooling neutron stars; e .g., 
lYakovlev fc PethicklEool iPaee. Geppert fc Weberll2006l) . 

Figure 2] demonstrates sensitivity of the results to the 
values of the thermal conductivity in the inner neutron star 
crust. It shows temperature profiles in the 1.4 Mq star of 
age 1000 years and the cooling curves of this star for one 
location of the heating layer (case I) and one heat intensity 
Ho = 3 x 10 20 erg cm -3 s -1 . The thick lines are the same 
as in Figure [5] They are calculated with our standard cool- 
ing co de assuming only electron thermal co nductivity in the 
crust iGnedin. Yakovlev fc Potekhinl200 However, in the 
inner crust thermal energy can also be transported by free 
neutrons and this transport can be very efficient, especially 
if neutrons are in a superfluid state. The effect may be sim- 
ilar to that in superfluid 4 He, where no temperature gradi- 
ents can be created in laboratory experiments because they 
are i mmediately smeared o ut by responding convective flows 
(e.g.. iTillev fc Tille^ll99ch . To simulate this effect we have 
artificially introduced the layer of high thermal conductiv- 
ity in the inner neutron star crust, in the density range from 
3 x 10 12 g cm" 3 to 10 14 g cm -3 . In this layer we enhanced 
the thermal conductivity by a factor of C — 10 4 or 10 s . The 
enhancement changes dramatically the temperature profiles 
in the inner crust, making it much cooler and almost isother- 
mal. We have also made one test run by reducing the thermal 
conductivity in the indicated layer by a factor of 10 4 . This 
corresponds to C = 10 -4 and mainly increases the temper- 
ature in the inner crust. However, all these strong changes 
of the temperature profiles in the inner crust have almost 
no effect on the surface temperature and the cooling curves. 
It is another manifestation of the thermal decoupling of the 
surface layers from the inner parts of the star. In addition, 
we have simulated the cooling of the star (the dotted curves) 
in the approximation of infinite thermal conductivity every- 
where in the star bulk (p > pb = 10 10 g cm -3 ). In this case 
the heat energy is instantly spread over the star bulk, which 
makes the stellar surface much cooler than in the case of 
finite conduction. 

Figure shows the integrated heating rate W°° , given 
by Eq. @ (three upper lines on each panel), and the photon 
thermal surface luminosity of the star (redshifted for 
a distant observer, three lower lines on each panel) as a 



function of parameters of the heating layers. The results are 
presented for the 1.4 Mq star of age 1000 years. 

In the left panel we select three locations of the heating 
layer (I, II, III) and vary the heat intensity Ho- One can 
clearly see that only the heating layer I can produce > 
3x 10 35 erg s -1 , typical for magnetars. Moreover, the surface 
luminosity increases with Ho much slower than the heating 
rate. For Ho > 10 20 erg cm -3 s -1 and the layers II and III, 
the luminosity is seen to saturate. This means that pumping 
additional heat energy into the heating layer does not affect 
L^°. The efficiency of converting the input heat into the 
surface radiation (L^° /W°°) is generally small. The highest 
efficiency is achieved if we heat the outer crust (the layer I) 
with low intensity. 

In the right panel of Figure|S]we present and W°° as 
a function of the maximum density p2 of the heating layer, 
for one heat intensity Ho = 10 20 erg cm -3 s -1 and three 
fixed minimum densities of this layer (pi = 3 x 10 10 , 10 12 , 
and 3 x 10 13 g cm -3 ). One can observe the saturation of 

with increasing p2. If p2 is large and the heating layer 
is extended into the core, the heat energy is huge. However, 
as long as pi is far from the surface, this huge energy is 
almost fully carried away by neutrinos and does not heat 
the surface. In this case the efficiency of heat conversion 
into the surface emission is very low. 

In Figure E| we compare qualitatively the calculated val- 
ues of with the thermal surface luminosities from the 
"magnetar box" (the lower shaded strip, estimated using 
accepted values of T a °° and the 1.4 Mq neutron star model). 
The heating should be sufficiently intense to raise to 
these magnetar values. It is the first requirement to explain 
the observations of the magnetars. The second requirement 
stems from the energy budget of an isolated neutron star. 
The heat energy is assumed to be taken from some source 
which pumps the energy W°° into the heating layer during 
magnetar's life (r~5x 10 4 years in Eq. Q). Naturally, the 
total available energy E to t is restricted. At least it should 
be much smaller than ~ 5 x 10 53 erg, the gravitational en- 
ergy of the neutron star. We assume that the maximum 
energy of the internal heating is -E ma x ~ 10 50 erg (which 
is the magnetic energy of the star with the magnetic field 
B — 3 x 10 1B G in the core). Then the maximum energy gen- 
eration rate is Wmax ~ Bmax/f ~ 3 x 10 37 erg s -1 , which is 
plotted by the upper horizontal solid line in Figure |S] The 
upper shaded space above this line is thus prohibited by the 
energy budget. 

A successful interpretation of magnetars as cooling neu- 
tron stars requires to be sufficiently high to reach the 
"magnetar box" and W°° to be sufficiently low to avoid 
the prohibited region. These conditions are fulfilled for the 
heating layer located in the outer stellar crust and the heat 
intensities Ho between 3 x 10 19 and 3 x 10 20 erg cm -3 s -1 
(higher value preferred). A typical efficiency of heat con- 
version into the surface emission under these conditions is 
L™/W°° ~ 10 -2 . 



4 DISCUSSION 

4.1 Testing the results 

We have performed a number of additional tests. In partic- 
ular, we have studied sensitivity of our results to the values 
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Figure 4. The effect of the thermal conductivity in the inner crust on the thermal evolution of the 1.4 Mq magnetar with the heating 
layer I and the heat intensity Hq = 3 x 10 20 erg cm -3 s — x . Left: The temperature profiles within the magnetar at t =1000 years. The 
hatched rectangles show the positions of the heating layer and the layer, where the thermal conductivity was modified. Right: The cooling 
curves. The thick solid lines are the same as in Figurc|5j Thinner long-dash, dot-and-dash, and short-dash lines are for the star with the 
thermal conductivity modified by a factor of C = 10 4 , 10 s , and 10 — 4 , respectively, in the density range from 3 X 10 12 to 10 14 g cm . 
The dotted line marked isothermal is for an infinite thermal conductivity in the star bulk (p > 10 10 g cm -3 ). 
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Figure 5. The total heat power W°° (higher curves) and the surface photon luminosity (lower curves) versus parameters of the 
heating layer compared to the values of from the "magnetar box" (the lower shaded strip) and to the values of W°° forbidden by 
energy budget (the upper shading) for the 1.4 Mq neutron star of age 1000 years. Left: Three fixed positions of the heating layer of 
variable heat intensity Hq. Right: Three fixed minimum densities pi of the heating layer, the fixed heat intensity Ho = 10 20 erg cm -3 s _1 
and variable maximum density p2- 



of the neutrino emissivity in the neutron star core and crust. 
We find that variations of the neutrino emissivity in the in- 
ner crust and the core of the star that undergoes intense 
heating can strongly affect the temperature profiles in the 
stellar interior but have almost no effect on the surface tem- 



perature. This is another example of the thermal decoupling 
of the outer crust and deep interiors. On the other hand, we 
find that is sensitive to the neutrino emission in the 
outer crust. This emission in a hot crust is mainly gener- 
ated by plasmon decay and electron-nucleus bremsstrahlung 
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mechanisms (see, e.g., lYakovlev et al-lhoOll) . If magnetars 
have very strong magnetic fields in their outer crusts, these 
fields can greatly modify the plasmon decay neutrino pro- 
cess. Such modifications have not been studied in detail but 
can be important for the magnetar physics. 

We have made some test runs taking into account the 
effects of superfluidity. In particular, we have included super- 
fluidity of free neutrons in the inner crust and the associated 
neutrino emission via Cooper pairing of neutrons. This neu- 
trino emission affects the temperature profiles in the inner 
crust but has almost no effect on the surface temperature. 

We have also varied the maximum density pb of the 
heat blanketing layer, shifting it from the present position 
p h = 10 10 g cm -3 (S ectional to p h = 4 x 1 11 g cm" 3 , as in 
our previous model jPotekhin et al J [20031. In addition, we 
have varied the magnetic field strength B in the blanketing 
envelope from B = 5 x 10 13 G to 5 x 10 15 G. The cooling 
curves and the internal temperature profiles are sensitive to 
these variations but do not violate our principal conclusions. 
It would be desirable to reconsider the cooling with a more 
careful treatment of heat transport and neutrino emission 
in a magnetized outer crust. 

4.2 The nature of internal heating 

The development of a specific theoretical model of the 
internal heating is outside the scope of this paper (ex- 
isting models are summa rize d, e.g., in review papers by 
IWoods fc Thompsonl200d and lHevl200r3) . Nevertheless, our 
results place stringent constraints on the possible models. 
The main requirement is that the total heat energy should 
be huge, -Etot ~ 10 49 — 10 50 erg, and should be released 
during 10 4 — 10 5 years in the outer neutron star crust. This 
does not necessarily require that the energy be stored in the 
outer crust but that it be effectively transformed into heat 
there. 

Our results agree with the widespread point of view that 
the magnetars can be powered neither by their rotation, nor 
by accretion, nor by thermal energy in a cooling star, nor by 
the strain energy accumulated in the crust. All these energy 
sources contain much less than 10 49 erg. 

Nevertheless, the required energy can be accumulated 
in the magnetic field if, for example, the star possesses the 
field B ~ (1 - 3) X 10 16 G in its core. The evolution of this 
magnetic field can be accompanied by a strong energy re- 
lease in the outer crust, where the electric conductivity is 
especially low and the field undergoes the strongest Ohmic 
dissipation. An actual structure of the magnetic field in the 
star can be complicated. In particular, the magnetic config- 
uration in the crust can strongly deviate from a magnetic 
dipole. 

One cannot exclude that thermal radiation of magne- 
tars is emitted from a smaller part of the neutron star 
surface (e.g., from hot spots near magnetic poles), im- 
plying lower total heat energies. Nevertheless, observable 
thermal X-r ay luminosities of mag netars ~ 10 34 - 10 36 
erg s" 1 (e.g..lMereghetti et alJliool iKaspi fc Gavriilll20o"l 
IWoods fc T hompson! 120061) are consistent with the interval 
of thermal luminosities in Figure calculated assuming the 
emission from the entire surface. Radiation from the sources 
included in our "magnetar box" usually demonstrates low 
pulsed fraction (< 20%) which also indicates that this radi- 
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ation can be emitted from the entire surface (although the 
pulsed fraction can be lowered by the gravitational lensing 
effect). 

Let us recall also alternative theories of magnetars. 
They suggest that the main e nergy release occurs in the 
magnetar's magnetosphere (e.g.. lBeloborodov fc Thompson! 
2006) and the radiation spectrum is formed there as a re- 
sult of comptonization and reprocession of the quasi-thermal 
spectrum. 



5 CONCLUSIONS 

We have modeled thermal states and thermal evolution of 
magnetars (SGRs and AXPs) with the aim to explain high 
surface temperatures of these neutron stars and their energy 
budget. Our main conclusions are as follows. 

(1) It is impossible to explain high thermal states of 
magnetars as cooling neutron stars without assuming that 
they undergo powerful heating. We have developed the idea 
that the heat source operates in the interior of magnetars. 

(2) The heat source can be located in a thin layer at 
densities p < 5 x 10 11 g cm -3 in the outer magnetar crust, 
with the heat intensity ranged from ~ 3 x 10 19 to 3 x 10 20 
erg cm -3 s _1 . The source cannot be located deeper in the 
interior because the heat energy would be radiated away 
by neutrinos; it would be unable to warm up the surface. 
This deeper heating is extremely inefficient and inconsistent 
with the energy budget of neutron stars. Pumping huge heat 
energy into the deeper layers would not increase the surface 
temperature. 

(3) Heating of the outer crust produces a strongly 
nonuniform temperature distribution within the star. The 
temperature in the heating layer exceeds 10 9 K, while the 
bottom of the crust and the stellar core remain much colder. 
The thermal state of the heat layer and of the surface is 
almost independent of physical parameters of deeper lay- 
ers (such as the equation of state, neutrino emission, heat 
transport, superfluidity of baryons), which means thermal 
decoupling of the outer crust from the inner layers. The to- 
tal energy released in the heat layer during magnetar's life 
(~ 10 4 - 10 5 years) cannot be lower than 10 49 - 10 50 erg; 
maximum 1% of this energy can be spent to heat the surface. 

The present calculations can be improved by a more 
careful treatment of heat conduction in a magnetized plasma 
of the outer crust, at p > 10 10 g cm' 3 , as well as by con- 
sidering different magnetic fields and the presence of light 
(accreted) elements in the surface layers of magnetars. We 
intend to study these effects in our next publication. 

The nature and the physical model of internal heating 
are still not clear; they should be elaborated in the future. In 
any case, one should bear in mind that the heat energy re- 
leased within the magnetars should be at least two orders of 
magnitude higher than the photon thermal energy emitted 
through their surface, and the energy release should take 
place in the outer stellar crust. These model-independent 
conclusions are consequences of the well-known principle 
that any hot, dense matter in stellar objects is a strong 
source of neutrino emission. 
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